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Abstract. We study the quantum mechanical Liouville model with attractive 
potential which is obtained by Hamiltonian symmetry reduction from the system 
of a free particle on SL(2, (R). The classical reduced system consists of a pair of 
^ ■ Liouville subsystems which are 'glued together' in such a way that the singularity 

of the Hamiltonian flow is regularized. It is shown that the quantum theory of this 
reduced system is labelled by an angle parameter 8 £ [0, 27r) characterizing the 
self-adjoint extensions of the Hamiltonian and hence the energy spectrum. There 
exists a probability flow between the two Liouville subsystems, demonstrating that 
the two subsystems are also 'connected' quantum mechanically, even though all the 
wave functions in the Hilbert space vanish at the junction. 
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1. Introduction 



The Liouville model has been under intense scrutiny in recent years, due to its relation to 
two dimensional quantum gravity, which is important in the theory of the world sheet of 
string theories. The model, which is solvable classically, has proved to be a rich source for 
developing techniques as well as for probing the universal features of quantum gravity in 
higher dimensions. However, the problem is that despite many fruitful achievements the 
model still resists a full understanding as a quantum theory. It is therefore heartening to 
observe that its toy model version obtained by ignoring the space dimension governed by 
the Hamiltonian, 

H(n,x) = 7T 2 + / ue 2x , (1.1) 

can be solved completely for \i > even quantum mechanically. The quantum mechanical 
Liouville model possesses a continuous energy spectrum and its eigenstates are given by 
modified Bessel functions [1] . A somewhat peculiar aspect of the quantum theory is that it 
has no vacuum state even though the energy is bounded from below, an aspect that stems 
from the repulsive exponential potential, which has no minimum. Now one might ask what 
happens if we replace the potential 'wall' with a 'well' by putting \i < 0. This attractive 
exponential potential will undoubtedly give rise to difficulty in quantization because the 
energy would then be unbounded from below without any bound states, indicating the 
quantum instability of the system. This is a reflection of the classical instability whereby 
the particle sinks indefinitely fast into the well, as the classical solution develops a sin- 
gularity and blows up at some finite time if ji < 0. This is perhaps the reason why the 
quantum theory of the Liouville model with attractive potential has not been considered 
seriously so far. 

Meanwhile, we have found in the study of W-algebras and the generalized KdV sys- 
tems that the Toda field theory, obtained by Hamiltonian symmetry reduction [2] from the 
Wess-Zumino-Novikov-Witten model based on a Lie group G, has a certain global struc- 
ture [3]. More precisely, the reduced theory is not merely a Toda theory but consists of 
a multiple of Toda theories as subsystems having both repulsive and attractive potentials 
in general. In particular, for G = SL(2, (R) we have a pair of Liouville models as subsys- 
tems in the reduced system, both of which have an either repulsive or attractive potential 
depending on the reduction performed. The interesting observation made there [3, 4] (see 
also [5]) was that, in the toy model version where the Toda field theories become the Toda 
lattices, the singular classical solutions that arise in the Toda lattices are regularized au- 
tomatically by the Hamiltonian reduction. An intuitive picture of the regularization may 
be gained by considering the simple case G = SL(2, (R) where we get two Liouville (toy) 
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models which are 'glued together' by identifying the limits x — > oo of the two models (see 
Fig. 1). The singular solution in one Liouville model is regularized by continuing it in time 
to the solution in the other Liouville model, causing the particle to oscillate between the 
two subsystems. This observation motivated us to investigate more closely the Liouville 
model with attractive potential, now regularized in the above sense, to see if any sensible 
quantization is possible. 

A first step in this direction was made in a paper by Fiilop [6], where, like the for- 
mer repulsive potential case, the theory is solved completely at the quantum level yielding 
Bessel functions as energy eigenstates. The salient result of [6] is that the spectrum is dis- 
crete, which is perhaps a natural consequence of the regularized classical solutions being 
oscillations, and that there are inequivalent quantizations characterized by certain param- 
eters specifying the self-adjoint quantum Hamiltonian and hence the spectrum. However, 
the argument in [6] appears to be unnecessarily complicated at a few crucial points, espe- 
cially when the self-adjoint Hamiltonians are constructed over the entire reduced system. 
The aim of this paper is to present a quantization approach which is much simpler and 
more direct in these points, and to furnish a complete version of the (regularized) quantum 
mechanical Liouville model with attractive potential. We shall find that, as in [6], there 
arise inequivalent quantizations but they can be characterized just by an angle parame- 
ter 9 G [0, 27r), and that the discrete energy spectrum obtained turns out to be different 
from that of [6]. We also see more naturally a probability flow between the two Liou- 
ville subsystems, a fact demonstrating that these subsystems are also connected quantum 
mechanically. 

The plan of the paper is as follows: To make the paper self-contained, in section 2 
we provide a necessary background for the classical reduced system. Then in section 3 
we present a quantum theory of the reduced system, which is a combined system of two 
Liouville models with attractive potential. The final section is devoted to discussion. We 
provide two appendices; Appendix A for a brief review of the general theory of self-adjoint 
extensions of symmetric operators, and Appendix B for a collection of formulae involving 
Bessel functions used in the text. 

2. Hamiltonian reduction and the global structure 

In order to set the scene for the system for which we discuss the quantization, we here 
recall the Hamiltonian symmetry reduction [2] which leads to the system of regularized 
Liouville models, together with the global description of the system developed recently [3] . 
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The reduction is the special case n = 2 of the Hamiltonian reduction that yields a 
multiple (2 n_1 ) of open, finite Toda lattices from the free particle system on the group 
G = SL(n, K). In the reduced system these Toda lattices, which have in general both 
repulsive and attractive potentials, are 'glued together' in such a way that no singularity 
arises. In this sense the reduction provides a natural means to regularize the singularities 
which exist in those Toda lattices that have attractive potentials. (This regularization is 
an example of a more general idea put forwarded originally in [5].) 

2.1. Classical Hamiltonian reduction 

The free particle on a semisimple Lie group G is described by the Hamiltonian system 
(M,Q,H) in the following way. The phase space M is the cotangent bundle of the group, 

M = T*G ~ G x Q = {(<7, J) \ g e G, J E G } , (2.1) 

where Q is the Lie algebra of G (in our case Q = sl(2, R)), which is identified with its dual 
Q* using the scalar product. The fundamental Poisson brackets are 

{9ij,9kl} = 0, {g %J MT a J)} = {T a g) l3 , {tr(T a J) , tr(T b J)} = tr([T a , T b ] J), (2.2) 

where {T a } is a basis of Q. For Q = s/(2, IR) we may take 

The Poisson brackets (2.2) derive from the symplectic form, 

O = dtr (Jdgg- 1 ) . (2.4) 



The Hamiltonian is given by 



which leads to the dynamics 



H = Uv(J 2 ) : (2.5) 



^ = {g,H} = Jg, ^ = {J,tf} = 0, (2.6) 

and hence yields the geodesic equation ^(^f <7 _1 ) = on the group manifold G. We note 
that J is the infinitesimal generator for the action of G on M defined by left translations, 
while the action of G defined by right translations is generated by J : M — > Q where 

J{g,J):=-g- 1 Jg. (2.7) 
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Let us now decompose the Lie algebra Q = sl(2, IR) into the subalgebras of strictly 
upper triangular, diagonal, and strictly lower triangular traceless matrices, that is, those 
subalgebras spanned by the basis (2.3), 

q = g+ + g + G-. (2.8) 

Then we consider symmetry reduction based on the subgroup N := N + x iV_ C G, where 
N± = exp(C?±), which acts on the phase space M according to 

(n+,n_) : (g, J) ^ (n+gnZ 1 ,n+Jn+ 1 ), V(n+,n_) e N, (<?, J) G M. (2.9) 

The symmetry reduction is performed by decomposing the generators, J = J + + Jq + J- 
and J = J_|_ + Jo + J-, according to (2.8) and then fixing the value of the momentum map 
$(g,J) := (J-,J+) as 

0>(<7, J) = (/_,-/+), (2.10) 

where 7_ := z/~T_ and I + := v + T + are nonvanishing constant matrices {v^ ^ 0) belonging 
to Q- and Q + , respectively. The reduced phase space is obtained as the factor space 

M red (/_, 1+) = M c (I_,I + )/N, where M C (I_,I + ) := (/_, -/+) . (2.11) 

In Dirac's terminology, this Hamiltonian reduction amounts to imposing the first class con- 
straints, J- = I- and J_|_ = defining M c C M, and getting the reduced phase space 
by fixing the gauge associated with the symmetry group iV generated by the constraints; 
hence (2.11). 

Now the Bruhat ( Gelfand-Naimark) decomposition for semisimple Lie groups [7] allows 
us to write G = SL(2, (R) as 

G = G e U G- e U Gi ow (disjoint union), (2.12) 

where 

G± e := ±N+AN_ with A := exp (<? ) • (2.13) 

The two 'cells', G e and G- e , are open submanifolds in G containing e and — e (e G G is the 
identity element), respectively (and their union is dense in G), while Gi ow is the union of 
'borders', i.e., certain lower dimensional submanifolds of G. Correspondingly, any element 
g G G e U G- e C SL{2, IR) admits the unique decomposition in the form, 

9 = ±^n.=±(l e l)(l J), ( 2 ,4) 



5 



where we put q = xT Q and a, c, x G R. The two cells G± e are in fact the open submanifolds 
of determinant one matrices with sign((7 2 2) = ±, whereas Gi ow consists of those matrices 
with (722 = 0. 

The Bruhat decomposition naturally induces the decomposition of the phase space 
M = T*G as M = M e U M_ e U M\ ow . Since this decomposition is invariant under the 
action of the symmetry group N, we have the corresponding decomposition of M c = 
Mg U Mt e U Mjq W as well, which in turn induces the decomposition of the reduced phase 
space, 

M red = M red y M red y M red (disjoint Union) . (2.15) 

We note that Mgf are open submanifolds in M red and Mj^ is a union of lower dimensional 
submanifolds. In other words, the Bruhat decomposition introduces the cell-structure in 
the reduced phase space. We now show that each of the subsystems associated with the 
two cells M^ e e d is a Liouville model. Indeed, since the submanifolds M± e are 

M c ±e = {(g,J)\g= ±n+e%_, n± G N±, q G Q , J_ = /_, (s" 1 J^) + = /+ } , (2.16) 
we see that M^r = M± e /(N + x iV_) are given by the local gauge section, 

M£ d = { (±e«, J) | q G So, J- = (e- fl Je«) + = J+ } . (2.17) 
The condition in (2.17) can easily be solved for J: 

J = J±e(q,p) := I-+p + e q I + e~ q , (2.18) 
where q, p G Go- Thus we may write 

Mg = { (±e«, J± e (q,p)) | ( ?J p) G Go x So } , (2.19) 

or simply M±£ ~ So x £o = R 2 - By evaluating the symplectic form (2.4) on the reduced 
phase space (2.19), we find 

O r ± ed = dtr (J ±e d(e q ) e~ q ) = dtr (pdq) = 2dndx, (2.20) 

where we put p = ttTq with tv G (R. Similarly, from (2.5) the reduced Hamiltonian turns 
out to be 

Hf d e {J±e) = itr (4 e ) = ^tr(p 2 ) +^tr(T_e^T + e^) = tt 2 (2.21) 
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where n = v + v~ . Since the triples (M± ed , n^r d , H±£) are no other than that of the Liou- 
vile model, we conclude that the reduced system obtained by the Hamiltonian reduction 
possesses two, identical Liouville models as subsystems. 

2.2. Global structure of the reduced system 

The Hamiltonian flow on the manifold M^ ed is governed by the reduced Hamiltonian (2.21), 
which yields the equation of motion, 

^ + t*e** = 0. (2.22) 

The flow is incomplete (singular) if n < 0, as is intuitively clear from the rapidly decreasing 
potential well in which the particle sinks. For instance, the solution, x(t) = — ln(cost) 
which satisfies the initial condition x(0) = and 4|(0) = and corresponds to \x = — 1, 
blows up at f = |. But since there is no singularity in the full reduced system, the 
incompleteness of the Hamiltonian flow that arises in the reduced system when \x < is 
just a manifestation of the fact that the particle may leave the submanifold M±°- C M red 
at finite time. More concretely, the trajectory of the free particle on G determined by 
an initial value (±e 9 , J± e ) at t = as g(t) = ±e tJ±e e q may leave the open submanifold 
G± e , because the flow of the reduced system is obtained by projecting the original flow 
on M c C M to M red . Thus the singularity occurs when (722 vanishes, which corresponds 
to q (or x) reaching infinity. In this respect, one can say that the embedding of Mj_ ed 
into M red provides a regularization of the Liouville model with \x < where the singular 
(blowing up) trajectories are glued together smoothly at x infinity. If, on the other hand, 
\i > 0, then the Hamiltonian flow is complete, and the two Liouville models are completely 
disconnected from each other. Hereafter we confine ourselves to the case \x < where the 
two Liouville models are connected, and put \i = — 1 for simplicity since \i can be freely 
rescaled by shifting x in (2.21). We also choose u + = —v~ = 1 for definiteness. 

So far, we have analyzed the structure of the reduced system only locally, using the 
gauge fixing to identify the reduced system (M red , O red , H ved ) as one containing a pair of 
Liouville models glued together along lower dimensional sub manifolds. To furnish a tool 
to gain information on the global structure, we wish to have a global cross section (gauge 
fixing) of the gauge orbits in M c . Such a cross section is furnished by the Drinfeld-Sokolov 
gauge, which is used in the context of generalized KdV equations [8] and W-algebras (see, 
e.g., [2]). In our context we need to use it doubly for J and J, and for this reason we call 
our global gauge fixing 'double DS gauge' here. 
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A double DS gauge is defined by requiring that the two generators J and J of the 
symmetry be of the dual form, 

J(t0 = (J «*), j (t0 = (£ J). ,, 23 , 

Note that by definition J and J are not quite independent (see (2.7)). It is readily seen 
that the condition for J in (2.23) is fulfilled if J is of the form in (2.23) and 

sm=(t :::)• <™> 

These parameters (112, u%, U4) G fR 3 are subject to the condition, 

F(u) := detg(u) = uj- u 2 ul = 1, (2.25) 

which defines the hypersurface S in K 3 as a model of M red . It is straightforward to check 
that dF(u)\F( u )=i 0, which implies that (2.25) gives a regular hypersurface diffeomorphic 
to M red . Regarding the ii; as gauge invariant functions on the constrained manifold M c , 
we find the Poisson brackets 

u 2 

{u 4 ,u 3 } = y, {u 4 ,u 2 } = us, {u 3 , u 2 } = u 4 u 2 . (2.26) 

The dynamics of the reduced system is determined by the Hamiltonian 

H red (u) = \trJ 2 {u) = u 2 . (2.27) 

The relationship to the local description of the reduced system given earlier is established 
by noting that, since #22 (w) = — U4, the cells M e red and Mi e e d are represented by the domains 
114 < and 114 > 0, respectively. 

Having obtained a global picture of the reduced system, we now see how the singularity 
of the Liouville model gets regularized when the two models are glued together. Consider 
the classical solution with a constant energy E, whose trajectory is the curve obtained by 
intersecting the hypersurface S in (2.25) with U2 = E. The curve is a hyperbola or ellipse 
depending on the sign of the energy E (see Fig. 2), 

u\ - Eu\ = 1. (2.28) 

It is now clear that, for E < 0, the motion of the particle is periodic, which implies that 
the particle does pass the border it 4 = along the curve, and when it does so it gives rise 
to a singularity in the solution when viewed as a local subsystem, although the solution is 
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perfectly regular when viewed as a global system. It is worth noting that the hypersurface 
S is not simply connected. In fact, the loop given by, say, the ellipse in (2.28) for some 
E < cannot be contracted to a point on the surface S. 



3. Quantization of the reduced system 



We have learned that the classical reduced system consists of two Liouville subsystems 
glued together between which the point particle oscillates if the energy is negative. In this 
section we wish to define a quantum mechanics of the reduced system and examine if this 
global feature appears at the quantum level as well. In quantum mechanics observables 
are represented by self-adjoint operators. Here we consider the self-adjoint Hamiltonian 
of our reduced system, which is perhaps the most fundamental observable, as a crucial 
ingredient to set up the quantum theory. The standard procedure [9] for finding self-adjoint 
Hamiltonians is to choose a suitable domain where the quantum Hamiltonian, now takes 
some operator form, becomes a symmetric operator, and then find an extended domain 
where it becomes a self-adjoint operator. (For a brief review of the general procedure, 
see Appendix A.) After this procedure, we shall find that there exists a probability flow 
between the two subsystems even though all the wave functions vanish at the junction 
of the two Liouville subsystems, a fact that illustrates that the two subsystems are also 
connected at the quantum level. 

3.1. Hamiltonian as a symmetric operator 

The basic problem for quantizing our reduced system is that the reduced phase space M red 
is not quite a cotangent bundle of some configuration space, but a system of two cotangent 
bundles nontrivially combined. To take this feature into account, we wish to formulate the 
quantum theory dealing with the two subsystems simultaneously, in such a way that the 
classical connectedness of the two subsystems will also be realized quantum mechanically. 

In section 2 we have seen that there exists a global description of the reduced system 
using the hypersurface S as a model for the reduced phase space. Among the three 
parameters used we choose the variable 



and regard it as 'coordinate' of the particle. The variable is chosen on the ground that the 
component #22 is gauge invariant under the symmetry action (2.9) and hence from (2.14) 
we have the direct relation Q = ±e~ x with the local Liouville coordinate x. The variable Q 
is convenient for our purpose since Q > corresponds to the subsystem (M e red , 0g 6d , Hl ed ) 




(3.1) 
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and Q < to (Mf_ e e d , f^ d , H r ^ ), respectively. The canonical momentum conjugate to Q 
is then given by P := —2^3/114 which satisfies 

{Q,P} = 1, (3.2) 

under the reduced Poisson brackets (2.26). In terms of these variables the classical reduced 
Hamiltonian (2.27) reads 

H =\Q 2p2 -^2- (3-3) 

Upon the identification P = 2e x rc the Hamiltonian (3.3) agrees with the local expression 
(2.21). Now the trouble is that the Hamiltonian is ill-defined at Q = 0, i.e., at the 
junction between the two subsystems. We exclude this point from the domain where Q 
is defined: Q G IR* := IR\{0} = [R + U [R~. As we shall see soon, despite this exclusion 
and the apparent trivialization of the reduced system into two decoupled subsystems, it 
is possible to construct a quantum theory such that the two Liouville subsystems are 
connected nonetheless. 

In quantization we elevate these canonical variables to linear operators on a Hilbert 
space with the Poisson bracket (3.2) replaced by the commutation relation, 

[Q,P]=i. (3.4) 

Working with the coordinate representation, we define the Hilbert space by the space of 
square integrable functions, 

H:={<f>\\\<f>\\<oo}, (3.5) 
where \\<fi\\ = a/ (</>, </>) is the norm of the state 0. We furnish the innerproduct by 

(0,V) = (0,V')- + (0,V') + , (3-6) 

where 

(0, v>>- := / ° H wrmi (0, := ^ HQrm)- (3.7) 

The measure dQ/\Q\ used in the innerproduct (3.7) can be derived, for instance, by the 
path-integral reduction, where the original measure f| t 3 (t) for the free particle system on 
SL{2, (R) reduces to the form Y[ t dQ/\Q\(t) after we integrate out the momentum variables 
J taking into account the constraints and the gauge fixing conditions [10]. A more direct 
way to see this is to consider the phase space path-integral for the reduced system with the 
Hamiltonian (3.3) and then integrate on P to get the configuration space path-integral, 
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which results precisely in the measure dQ/\Q\. Note that the measure dQ/\Q\ is just the 
standard Toda measure dx in terms of the local coordinate x, as expected. Note also that, 
because of the measure, all the wave functions <p{Q) G H must vanish at the junction of 
the two subsystems: 4>{Q) — > as Q — > ±0. 

In order to find a self-adjoint Hamiltonian operator, let us consider the differential 
operator H of the form, 

AAA 1 

(3.J 



which is a naive choice for the operator which corresponds to the classical Hamiltonian 
(3.3). We can find a domain where the Hamiltonian operator is symmetric, 



D(H) := {iplip eH, Hip e H, lim Qip{Q) 
In fact, it can be readily confirmed that on D(H) we have 



0}. 



(3.9) 



(3.10) 



It is also easy to see that the domain of the adjoint operator H*, which as a differential 
operator takes the same form as H 1 is just D{H*) := {tp | ip e H, Hip e H}. This shows 
that the symmetric operator H is not self-adjoint; D(H*) D D{H). 

Before discussing the self-adjoint extensions of the symmetric operator, let us consider 
the eigenvalue problem of the differential operator H, 



Using the variable 

we find that eq.(3.11) becomes 



H<f> = E<f>. 



d 2 d 



where 



z 2 —r + z— + (z 2 - k 2 ) 
dz z dz 



k 2 = -4E. 



= 0, 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



The linearly independent eigenstates of eq.(3.11) are given by the Bessel functions Jk{z) 
and Yk(z) with indices k E (C. However, only those Bessel functions of the type Jk{z) with 
Re k > have finite norms and belong to the Hilbert space H. In particular, for E < 
there exists a unique k for which the corresponding eigenstate belongs to H while there 
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is no such eigenstate belonging to H for E > 0, and for this reason we henceforth restrict 
ourselves to the negative energy states. Note however that none of these negative energy 
eigenstates belongs to the domain D(H) where the Hamiltonian is symmetric, as can be 
seen from their asymptotic forms (see Appendix B), 

2, ^ 1 ^3 ,2 tt„ 1. 



QMq) Q Z + o 0F^cos(---(/c+-)), 



„ d t ,2, „ 1 i ,2 tt /7 1,. ^ 3 ' 15 ' 1 

-^^"(g-^+j))- 

Clearly, as Q tends to zero the first term QJk(2/Q) goes to zero while the second term 
Q-^Jk(2/Q) blows up to infinity, showing that Jk{z) & D(H) for Rek > 0. 

3.2. Self-adjoint extensions 

In order to find an extended domain where the symmetric operator H becomes a self- 
adjoint operator H* for which D(H*) = D(H*), we first examine, according to the general 
theory (see Appendix A), the deficiency indices (d+,d-) of the symmetric operator H. 
The index d+ (d-) is given by the dimension of the eigenspace of the adjoint operator H* 
with eigenvalue % (—i), 

d+ := dimKer(^* -i), d- := dimKer(if* + i). (3.16) 

We recall that the eigenfunctions of the differential operator H in (3.8) with eigenvalue 
i (— i) are given by the Bessel functions J±k (J±k*) with index k := 2exp( — ji). But 
since for these eigenfunctions to be in the Hilbert space H the real part of the indices ±ko 
(±/cq) must be positive, the only eigenstate allowed is Jk (Jk*), which actually belongs to 
the domain of the adjoint operator D(H*). We therefore observe that d_) = (1, 1), a 
result ensuring that H admits self-adjoint extensions. 

The general theory of self-adjoint extension then asserts that the extended domain 
where the symmetric operator becomes self-adjoint consists of three spaces, i.e., the domain 
of the symmetric operator D(H), the eigenspace of the adjoint operator H* with eigenvalue 
i, and the eigenspace of H* with eigenvalue — i, with the latter two spaces being related 
unitarily. Concretely, it is given by 

D(H*) = {^\^ = ( p + aCfco, e D(H), a G 1}, (3.17) 

where 

{h -ra +e Par ( ' 
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The point to be noted is that the domain is parametrized by the angle, 



#e[0,27r). (3.19) 

The angle parameter 9 therefore characterizes the quantum theory we construct, whose 
existence is just a reflection of the ambiguity in quantizing a classical system. 

Having achieved the self-adjoint extensions of the Hamiltonian, we now investigate 
which eigenstates Jk of the adjoint operator belong to D(H*) in (3.17) for a specific 9. A 
necessary and sufficient condition for Jk G D(H*) is [11] 

(H*J k , iP) = (J k , H*ip) Vip £ D(H*). (3.20) 

which takes the form 



-o 



+00 



QW(Ji^) -QW(Ji^) =0, (3.21) 

— oo +0 

where W((pi, fa) = ^xjfa — fa^jQ is the Wronskian. An equivalent condition is obtained 
if we replace tp with £k in (3.21), and by using the asymptotic forms (3.15) we arrive at 
the relation 

e 2m(k-Rek ) = cos | cosh (Trim fc ) -zsinj sinh (Trim fc ) ^ 
cos | cosh (Trim fco) + i sin | sinh (Trim ko) 

Solving this relation in favour of k, we get 



k = k(n, 9) := Reko tan 1 

TT 



Q 

tan - tanh (Trim ko) 



+ n, (3.23) 



where n are integers for which k > 0. We therefore see that the eigenstates allowed by 
the self-adjoint Hamiltonian labelled by 9 are characterized by the indices k(n,9) whose 
energy eigenvalues are discrete, and that from (3.14) the intervals of adjacent discrete 
energy levels depend on the angle parameter. Note that as we vary the parameter from 
9 = the spectrum changes accordingly and returns to the original spectrum only when 
9 approaches to 2tt. Since any two Bessel functions whose indices differ by an integer are 
orthogonal to each other with respect to the innerproduct on IR*, so are any of the two 
eigenstates in the domain D(H*), as required. 

3.3. Probability flow between the two Liouville subsystems 

We have seen in section 2 that the reduced classical system admits solutions Q(t) oscillating 
between the two subsystems given by Q > and Q < 0. We now analyze how this global 
aspect manifests itself in the quantum theory. More specifically, we are interested in the 
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question if there exists a probability flow between the two subsystems. We shall find that 
the answer is positive, signaling the fact that the reduced system is a connected system 
also quantum mechanically. 

For this, consider the state ip(t) given by a linear combination of two energy eigenstates 

m = ciJk 1 (^)e- iEk ^ + C2Jk 2 (^)e- iE ^ t , (3.24) 

where J kl , J k2 G D(H*) for some fixed 9 and E ki = —kf/4 for i = 1, 2. We take the two 
indices k± and k 2 which are different by some odd integer, 

k 1 -k 2 = 2n+l, n G Z, (3.25) 

and choose the constants ci, c 2 G (C so that the state be normalized H^H = 1 at t = 0. Then 
the orthogonality condition (J kl , Jk 2 ) = implies that the norm, i.e., the total probability 
on the full line R*, remains constant, 

However, the probability on a half line, say K + , does not remain constant. Indeed, a similar 
computation for the half line reveals that 

j t (iP, V)+ = i4c 2 (J kl , J k2 ) + (E kl - E^e-W*'-^ + c.c. ^ (3.27) 

because the two Bessel functions with (3.25) are not orthogonal to each other (J kl , Jk 2 )+ 7^ 
on the half line (R + . This shows that there exists a probability flow between the two 
subsystems where Q G R and Q 6 R~, even though the wavefunctions vanish at the 
junction. The reason why such a flow can exist is that the Hamiltonian operator is not 
self-adjoint with respect to the innerproduct on the half line [R + (or (R ), although it 
becomes so if we cut the domain 'in half so that it consists only of those eigenstates with 
indices given by even (or odd) integers n in (3.23). 

4. Discussion 

We have seen in this paper that the quantum mechanical Liouville model with attractive 
potential obtained by the Hamiltonian reduction (which regularizes the Liouville model 
classically) can be solved completely. Although the energy spectrum is unbounded from 
below, the fact that only discrete levels are allowed suggests that the system is 'quasi- 
stable' at the quantum level. The connectedness of the two subsystems can be observed 
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by a probability flow, which we have shown to exist. The quantum theory is labelled by 
the angle parameter 9 which arises in constructing self-adjoint Hamiltonian operators. In 
this respect, it is worth noting that the reduced phase space is topologically IR 2 \{0}, that 
is, the two dimensional Euclidean plane with a hole, which is homeomorphic to T*S 1 = 
S 1 x IR. The appearance of the angle parameter may perhaps be understood as a common 
phenomenon observed in quantizing on a configuration space with a hole, as in the case 
of the quantum theory on S 1 or of the Yang-Mills theory, whose quantization yields the 
6>-vacua by an analogous mechanism. 

Our procedure of quantization is similar to that of Fiilop [6] but differs in some impor- 
tant points. First, unlike our measure dQ/\Q\ for the innerproduct (3.7), the measure used 
for the innerproduct in [6] is \Q\dQ, which is obtained from the Haar measure of the group 
SL{2, IR) by eliminating the degrees of freedom that correspond to the symmetry. However, 
in Hamiltonian reduction one has to use the phase space volume element on T*SL(2, IR) 
to derive the correct reduced measure [10], which is the measure dQ/\Q\ we used. This 
causes a certain alteration in the energy spectrum. Second, and more importantly, the 
self-adjoint extensions of the Hamiltonian operator are achieved in [6] by considering a 
domain such that the eigenstates on the full line IR* are formed out of linear combinations 
of two eigenstates, each defined on the half lines IR + and IR~ , respectively. Since one can 
take distinct angle parameters to specify the self-adjoint extensions on the two half lines, 
one needs two angle parameters in general to specify the self-adjoint extensions on the full 
line IR* (plus an extra parameter to render the eigenstates mutually orthogonal) . This we 
find an unnecessary complication, given that the self-adjoint extensions can be achieved on 
the full line without referring to those on the half lines. Our simpler quantization yields 
the energy spectrum given by a single class of discrete levels with indices (3.23) specified 
by the angle 9, a result which we feel is natural to associate with the classical periodic 
motions on the smooth phase space. In contrast, the spectrum in [6] consists of two classes 
of discrete levels similar to (3.23) but with even integers n G 2Z. 

Finally, we wish to stress that the quantization discussed in this paper is not the 
unique one available to the Liouville system. Indeed, from the way the classical system is 
defined, it is perhaps more natural to quantize first the system of a free particle on G and 
then carry out quantum Hamiltonian reduction. This will provide a way to confirm what 
we have learned in the quantum mechanical Liouville model given in this paper. 
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Appendix A. Symmetric operator and self-adjoint extension 

Here we briefly summarize the general theory of self-adjoint extensions of symmetric op- 
erators [9]. Let A be a linear operator on a dense domain D(A) in a Hilbert space TC. 
Consider if; G Ti for which 

(1>M) = WA) V^D(i). (A.l) 

is satisfied for some ip' G H. Define the domain D(A*) by the set consisting of those ip. 
It is then easily confirmed that ip > ip' gives a linear map. The adjoint operator A* is 
defined by this linear map A*: ip' = A*ip. If the operator A fulfills the condition 

D(A*) D D(A), A* = A on D(A), (A.2) 

then A is a symmetric operator. A symmetric operator is called a self-adjoint operator 
when the two domains coincide, D(A*) = D(A). If A, B are two operators such that 

D(B) D D(A), B = A on D(A), (A.3) 

then the operator B is an extension of the operator A. A symmetric operator A can be 
extended to be a self-adjoint operator A* if there exists a domain D(A Sf ) such that 

D(A*) D D(Al) = D(A*) D D(A). (A.4) 

When a certain condition is fulfilled (which we discuss shortly), such self-adjoint ex- 
tensions of a symmetric operator A are possible in the following way. We begin by decom- 
posing the domain D(A*) of the adjoint operator A* as 

D(A*) =D(A) +K-(A*) +K+(A*), where K±{A*) := Ker (A* ± i). (A.5) 

To see that this decomposition is possible, we first decompose any state ip G D(A*) as 
ip = ( + a where ( G D(A) and a G" D(A). Applying (A* — i) to ip we find 

(A*-i)ip= {A - i)( + (A* - i)a. (A.6) 

Note that on account of the property ||(A ± z)^>|| 2 = ||A^|| 2 + \\(f)\\ 2 for V</> G D(A) the 
spaces (A ± i)D(A) are closed subspaces in Ti, and that they are orthogonal to K^(A*), 
respectively. Thus we can write (A* — i)a = (A — i)(3 — 2z^, where /5 G -D(A) and ^ G 
K+(A*). But since -2^ = (A* - i)^, we have 

(A*-z)(^-0-C) = O, (A.7) 
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where := ( + (3 G -D(A). It then follows that rf :=ip — <f> — £ G if_(A*), that is, any state 
V> G _D(A*) can be decomposed as 

V> = + £ + *7, 0GD(i), 77GK + (i*), (A.8) 

which proves our claim (A. 5). 

It can be shown that a necessary and sufficient condition for satisfying the relation 

(V,i» = (i*V» (A.9) 

for ip G D(A*) in the form (A.8) is ||£|| = \\rj\\, i.e., the latter two components must be 
related by a unitary transformation rj = U£. Such a unitary transformation exists if and 
only if the two deficiency indices (d + , d-) defined by 

d+ = dimK_(i*), d- = dimif + (i*), (A.10) 

are equal d+ = d-. Hence, if this is the case, the symmetric operator A can be extended 
to a self-adjoint operator A* with the domain 

£>(i*) = {V>|V = + £ + t>£, 0e£>(i), £GK_(i*), U£eK + (A*)}. (A.n) 



Appendix B. Orthogonality and the asymptotic forms of Bessel functions 

Bessel functions Jk(z) are defined by the series, 

^ (-l) m s z xk+2m 

Jk{z) = -rFTTT —— TT U ' arg |^| < tt, fcet. (B.l) 

^— ' m! T (k + m + 1 \2/ 



m=0 



The other type of Bessel functions are Yk(z) := (cos(kiv)Jk(z) — J-k(z))/ sin(/c7r). The 
Bessel functions (B.l) have the following asymptotic forms for z approaching infinity and 
zero: 

/ 2 TV 

Jk(z) 2 ~ +00 \ — cos[z- -(fc + 1)], 

V nz 2 , B 2) 

T(\ ~ 1 (z\ k 

Jk[Z) z - +0 F(k + 1) \2J ' 

We define the Bessel functions on IR~ by analytically continuing them from the upper half 
plane Imz > 0, that is, we give the values on IR~ by rotating the functions on IR + by n, 

Jk(z) = Jk^K-z)) = e k ™J k (-z), wgz = jr. (B.3) 
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Note that JjJ(z) = J k *(z) for z£fl? + whereas J£(jz) = e- 2ik * n J k *(z) for 2 e IR". 

Due to the formula (B.3), integrals of the Bessel functions on the full line K\{0} read 

(J k ,J l ) = (l + e- i ^- l ^)(J k ,J l ) + . (B.4) 

The integrals ( J k , Ji) + on the half line K + can be evaluated by means of Lommel's integral 
[12], 



dz 



J k {z)Ji(z) 



1 



k 2 -l 2 



z(j k (z)^-Ji(z) - Ji{z)^-J k {z)) 



d 



dz 



dz 



(B.5) 



If Re (k + I) > then the above integral converges, and from Q = 2/z we obtain 



(B.6) 



b— >+oo 



where we used the asymptotic forms (B.2). If k, I e IR and /c — / = n G Z then the integral 
on the full line vanishes, (J/-, J{) = 0. This occurs for two different reasons depending on 
whether n is even or odd; for even n, the integral on the half line vanishes (J k , Ji)+ = 0, 
whereas for odd n the integral is not zero but the factor (1 + e~ l ^ k *~ l ^>) vanishes. 
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Figure Captions 



Figure 1. A schematic picture of the potential in the regularized Liouville system. 

Figure 2. Two typical trajectories on the reduced phase space. The ellipse (E < 0) 
corresponds to a periodic motion whereas the hyperbola (E > 0) corresponds to a motion 
passing once for all from one cell to another. 
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u 2 = - E ik 



